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Correlations among hadrons with the same electric charge produced in TA decays 
are studied using the high statistics data collected from 1991 through 1995 with 
the OPAL detector at LEP. Normalized factorial cumulants up to fourth order are 
used to measure genuine particle correlations as a function of the size of phase 
space domains in rapidity, azimuthal angle and transverse momentum. Some of 
the recently proposed algorithms to simulate Bose-Einstein effects, implemented in 
the Monte Carlo model Pythia, reproduce reasonably well the measured second- 
and higher-order correlations between particles with the same charge as well as 
those in all-charge particle multiplets. 


1 Introduction 

Correlations in momentum space between hadrons produced in high energy interactions have been extensively 
studied over many decades in different contexts^ Being a measure of event-to-event fluctuations of the number 
of hadrons in a phase space domain of size A, correlations provide detailed information on the hadronisation 
dynamics, complementary to thaC derived from inclusive single-particle distributions and global event-shape 
characteristics. The suggestion inH that multiparticle dynamics might possess (multi-)fractal properties or be 
“intermittent”, emphasized the importance of studying correlations as a function of the size of domains in 
momentum space. A key ingredient for such studies is the normalized factorial moment and factorial cumulant 
technique. Unlike factorial moments, cumulants of order q are a direct measure_pf the stochastic interdependence 
among groups of exactly q particles emitted in the same phase space cell.0“U Therefore, they are well suited 
for the study of true or “genuine” correlations between hadrons and are particularly sensitive to Bose-Einstein 
correlations. 

Two-narticle Bose-Einstein correlations (BEC) have beeiu observed in a wide range of multihadronic 
processesJj Such correlations were extensively studi^|at LEP.0 “eI Evidence for BEC among groups of more 
than two identical particles has also been reported.^^Uldl The subject has acuaLMed particular importance in 
connection with high-precision measurements of the VU-boson mass at LEP-ILllSlld For these, better knowledge 
of correlations in general is needed, as well as realistic Monte Carlo modelling of BEC. 

The high statistics OPAL data collected at and near the Z° centre-of-mass energy have been used to measure 
cumulants for multiplets of particles with the same charge, hereafter referred to as “like-sign cumulants”. They 
are compared to “all-charge” cumulants, corresponding to multiplets comprising particles of any (positive pc 
negative) charge. The role of Bose-Einstein-type effects is studied, using recently pmposed BEC algorithmalJ 
implemented in the Monte Carlo event generator Pythia for e+e” annihilation.td Proceeding beyond the 
usual analyses of two-particle correlations, we show that, at least within the framework of this model, a good 
description can be achieved of the factorial cumulants up to fourth order in one-, two- and three-dimensional 
phase space domains. 

2 The method 

To measure genuine multiparticle correlations in multi-dimensional phase space cells, w^use the technique of 
normalized factorial cumulant moments, AA, (“cumulants” for brevity) as proposed in.B The cumulants are 
computed as in a previous OPAL analysis.EA A D-dimensional region of phase space is partitioned into 
cells of equal size A. From the number of particles counted in each cell, rim {m = 1,. ■ ■ ,M^), event-averaged 
unnormalized factorial moments, (rim), and unnormalized cumulants, kq"^\ are derived, using the relations 
given e.g. in3. For g = 2, 3,4, one has 
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Here, (Til'll) = (n(n — 1 )... (n — g + 1)) and the brackets (•) indicate that the average over all events is taken. 
Normalized cumulants are calculated using the expression 

= (4) 

As proposed ir0, this form is used to correct for statistical bias and non-uniformity of the single-particle spectra. 
Here, is the number of particles in the mth cell summed over all N events in the sample, j■ 

The horizontal bar indicates averaging over the cells in each event, (1/M'°) 

Here, data are presented for “all-charge” and for “like-sign” multiplets. For the former, the cell-counts 
are determined using all charged particles in an event, irrespective of their charge. For the latter, the number 
of positive particles and the number of negative particles in a cell are counted separately. The corresponding 
cumulants are then averaged to obtain those for like-sign multiplets. 

3 Experimental details 

The analysis uses a sample of approximately d-Lx 10® hadronic Z® decays collected from 1991 through 1995. The 
OPAL detector has been described in detail in.lla The results presented are mainly based on the information from 
the central tracking chambers. The event selection criteria are based on the multihadronic selection algorithms 
described in.ll3 Multihadron events were selected with at least 5 good tracks, a momentum imbalance (the 
magnitude of the vector sum of the momenta of all charged particles) of less than 0.4 -y/i and the sum of the 
energies of all tracks (assumed to be pions) greater than 0.2 -^/i. In addition, the polar angle of the event 
sphericity axis, calculated using tracks and clusters had to satisfy | cos0gpjj| < 0.7 in order to accept only events 
well contained in the detector. A total of about 2.3x10® events were finally selected for further analysis. 

The cumulant analysis is performed in the kinematic variables rapidity, y, azimuthal angle, 4>, and the 
transverse momentum variable, InpT, all calculated with respect to the sphericity axis. 

• Rapidity is defined as y = 0.51n[(E + p^\)/{E — Py)], with E and py the energy (assuming the pion mass) 
and longitudinal momentum of the particle, respectively. Only particles within the central rapidity region 
—2.0 ^ p ^ 2.0 were retained. 

• In transverse momentum subspace, the logarithm of px is used to eliminate as much as possible the strong 
dependence of the cumulants on cell-size arising from the nearly exponential shape of the p|.-distribution. 
Only particles within the range —2.4 < In(pT) ^ 0.7 (pr in GeV/c) were used. 

• The azimuthal angle ( 0 ^ < 27r), is calculated with respect to the eigenvector of the momentum 

tensor having the smallest eigenvalue in the plane perpendicular to the sphericity axis. 

4 Results 

The fully corrected normalized cumulants Kq {q = 2,3,4) for all-charge and like-sign particle multiplets, calcu¬ 
lated in two-dimensional y x $ (2D) and three-dimensional y x $ x Inpr (3D) phase space cells, are displayed 
in Fig. 1^. It is seen that positive genuine correlations among groups of two, three and four particles are present: 
Kq > 0. Cumulants in 2D and 3D continue to increase towards small phase space cells. Moreover, the 2D 
and 3D cumulants are of similar magnitude at fixed M, indicating that the contribution from correlations in 
transverse momentum is small. The like-sign cumulants increase faster and approach the all-charge ones at 
large M. As the cell-size becomes smaller, the rise of all-charge correlations is increasingly driven by that of 
like-sign multiplets. 

The cumulant data have been compared with predictions of the Pythia Monte Carlo event generator 
(version 6.158) without and with Bose-Einstein effects. The model parameters, not related to BEC, were set 
values obtained from a previous tune to OPAL data on event-shape and single-particle inclusive distributions.Ej 
In this tuning, BE-effects were not included. 

To assess the importance of BE-type short-range correlations between identicaiuparticles, and their influ¬ 
ence on all-charge cumulants, we concentrate on the algorithm BE32, described incj, using parameter values 
PARJ(93) = 0.26 GeV (i? = 0.76 fm) and PARJ(92) = A = 1.5. Non-BEC related model-parameters were set at 
the following values: PARJ(21)=0.4 GeV, PARJ(42) =0.52 GeV"^, PARJ(81)=0.25 GeV, PARJ(82)=1.9 GeV. 
We find that calculations with PARJ (93) in the range 0.2 — 0.3 GeV, and the corresponding A in the range 
1.7 — 1.3, still provide an acceptable description of the second-order like-sign cumulants. 

The dashed lines in Fig. [| are Pythia predictions for like-sign multiplets for the model without BEC. 
Model and data agree for small M (large phase space domains), indicating that the multiplicity distribution in 
those regions is well modelled. However, for larger M, the predicted cumulants are too small. 
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Figure 1. The cumulants Kq in two-dimensional Ay X A<I> (2D) and three-dimensional Ay X A# X Alnpy (3D) domains for 
all charged hadrons (solid symbols) and for multiplets of like-sign particles (open symbols), versus M. Where two error-bars are 
shown, inner ones are statistical, and outer ones are statistical and systematic errors added in quadrature. Tihn lines connect Monte 
Carlo predictions from Pythia without BEC (dashed) and with BEC (full) simulated with algorithm BE32lifl (see text). 


The solid curves in Figs. show predictions for like-sign multiplets using the BE 32 algorithm. Inclusion 
of BEC leads to a very significant improvement of the data description. Also two-particle and higher order 
correlations in ID rapidity space are well accounted for (not shown). The predicted 2D and 3D cumulants agree 
well with the data. 

The ID, 2D and 3D cumulants for particle pairs with the same charge are displayed in Fig. Since BEC 
occur only when two identical mesons are close-by in all three phase space dimensions, projection onto lower¬ 
dimensional subspaces, such as rapidity and azimuthal angle, leads to considerable weakening of the effect. 
Nevertheless, the high precision of the data in Fig. || allows to demonstrate clear sensitivity to the presence or 
absence of BEC in the model. 

Whereas the BE-algorithm used implements pair-wise BEC only, it is noteworthy (see Fig. ^ that the 
procedure also induces like-sign higher-order correlations of approximately correct magnitude. This seems to 
indicate that high-order cumulants are, to a large extent, determined by the second-order one. 

To assess the sensitivity of the cumulants to variations in the BEC algorithms available in Pythia, we have 
further considered the algorithms BEa and BEoiiil Using the same parameter values as for BE 32 , we observe 
that BEa slightly underestimates K 2 {y) and overestimates K 2 {^) for like-sign pairs (Fig. whereas the results 
coincide with those from BE 32 in 2D and 3D. For all-charge multiplets (Fig. H), the predicted cumulants generally 
fall below those for BE 32 , except for and in 2D and 3D, where the differences are small. The differences 
with respect to BE 32 are related to the different pair-correlation functions used in the algorithms. Although a 
different choice of the parameters R and A may improve the agreement with the data, we have not attempted 
such fine-tuning. i—. 

We also considered the predictions based on the algorithm BEJIa (dash-dotted curves in the figures) for 
the same parameter values as quoted above. For like-sign pairs (Fig. ||), K 2 {y) and especially K 2 {^) are 
overestimated. In contrast, all-charge higher-order cumulants differ little from those obtained with BE 32 . 












Figure 2. The cumulants K 2 for like-sign pairs in one-dimensional domains of rapidity (y) and azimuthal angle (^), and in two- 
dimensional Ay X A<I> (2D) and three-dimensional Ay X A<I> X A InpT (3D) domains versus M. The error-bars show statistical and 
systematic errors added ini quadrature. The lines connect Monte Carlo predictions from Pythia, without BEC and with various 
Bose-Einstein algorithms^ (see text). 


To summarize, a comparison with Pythia predictions shows that short-range correlations of the BE-type 
are needed, at least in this model, to reproduce the magnitude and the A-dependence of the cumulants for like- 
sign multiplets. This further leads to a much improved description of the cumulants for all-charge multiplets. 
Since Bose-Einstein correlations are a well-established phenomenon in multiparticle production, it is likely that 
the above conclusion has wider validity than the model from which it was derived. 

The success of the Pythia model with BEC in predicting both the magnitude and domain-size dependence 
of cumulants, has led us to consider the inter-dependence of these quantities. Figure. || shows K 3 and K 4 
in 2D and 3D, as a function of K 2 . The 2D and 3D data for all-charge, as well as for like-sign multiplets 
follow approximately, within errors, the same functional dependence. The solid lines are a simple fit to the 
function XnKq = aq rq\YiK 2 . Figure ^ suggests that the cumulants of different orders obey simple ^ 
called “hierarchical” relations, analogous to the Ochs-Wosiek relation, hrst established for factorial moments 
Interestingly, all-charge as well as like-sign multiplets are seen to follow, within errors, the same functional 
dependence. 

Simple relations among the cumul^ts of different orders exist for certain probability distributions, such 
as the Negative Binomial distribution.Ell For this distribution, one has Kq = {q — {q = 3,4,...), 

showing that the cumulants are here solely determined by K 2 . This relation, shown in Fig. (dashed lines) 
does not describe the data, indicating that the multiplicity distribution of charged particles, and that of like-sign 
particles, deviates strongly from a Negative Binomial in small phase space domains. 

The Ochs-Wosiek type of relation exhibited by the data in Fig. ^ may explain why the BE algorithms in 
Pythia generate higher-order correlations of (approximately) the correct magnitude. 

5 Summary 

A comparative study of like-sign and all-charge genuine correlations between two and more hadrons produced 
in e+e" annihilation at the energy has been performed by OPAL using the high-statistics data on hadronic 
Z° decays recorded with the OPAL detector from 1991 through 1995. Normalized factorial cumulants were 
measured as a function of the domain size. A, in D-dimensional domains (D = 1,2,3) in rapidity, azimuthal 
angle and (the logarithm of) transverse momentum, dehned in the event sphericity frame. 
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Figure 3. The cumulants Kq in two-dimensional Ay X A<I> (2D) and three-dimensional Ay X A# X Alnpy (3D) domains for all 
charged hadrons versus M. Where two error-bars are shown, inner ones are statistical, and outer ones are statistical and systematic 
errors addedi in quadrature. The lines connect Monte Carlo predictions from Pythia, without BEC and with various Bose-Einstein 
algorithmsEj (see text). 


Both all-charge and like-sign multiplets show strong positive genuine correlations up to fourth order. The 
2D and 3D cumulants and K^, considered as a function of K 2 , follow approximately a linear relation of 
the Ochs-Wosiek type: \nKg ~ lniir 2 , independent of D and the same for all-charge and for like-sign particle 
groups. 

The Pythia model describes well dynamical fluctuations in large phase space domains. However, to achieve 
a more satisfactory data description, short-range correlations of the Bose-Einstein type between identical par¬ 
ticles need to be included. 

The Bose-Einstein model BE 32 in Pythia is able to simultaneously account for the magnitude and A- 
dependence of like-sign as well as of all-charge cumulants. The models BEq and BEa, when using the same 
parameters as for BE 32 , show reasonable agreement with the data. Although the algorithms implement pair-wise 
BEC only, surprisingly good agreement with the measured third- and fourth-order cumulants is observed. 
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